Thermodynamic properties of tetrameric bond-alternating spin chains 
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Thermodynamic properties of a tetrameric bond-alternating Heisenberg spin chain with 
ferromagnetic-ferromagnetic-antiferromagnetic-antiferromagnetic exchange interactions are studied 
using the transfer-matrix renormalization group and compared to experimental measurements. The 
temperature dependence of the uniform susceptibility exhibits typical ferrimagnetic features. Both 
the uniform and staggered magnetic susceptibilities diverge in the limit T 0, indicating that the 
ground state has both ferromagnetic and antiferromagnetic long-range orders. A double-peak struc- 
ture appears in the temperature dependence of the specific heat. Our numerical calculation gives 
a good account for the temperature and field dependence of the susceptibility, the magnetization, 
and the specific heat for Cu(3-Clpy)2(N3)2 (3-CIpy=3-Chloroyridine). 

PACS numbers: 75.10.Jm, 05.50.+q 



I. INTRODUCTION 

Physical properties of one-dimensional quantum fer- 
rimagnets have attracted much interest in recent years. 
In particular, physical properties of the so-called Lieb- 
Mattisi^type ferrimagnets have been extensively inves- 
tigated both theoretically and experimentallyiSiSi^ In 
these materials, there exists strong competition be- 
tween ferromagnetic and antiferromagnetic fluctuations. 
This competition leads to a ferrimagnetic-ordered ground 
state with coexisting gapless ferromagnetic and gapped 
antiferromagnetic excitations, a minimum in the temper- 
ature dependence of the product of the susceptibility and 
temperature, a double-peak structure in the specific heat, 
and a wide variety of other physical phenomena. 

A Lieb-Mattis-type ferrimagnet can be constructed 
from monospin chains with polymerized exchange inter- 
actions. It can be also constructed from an alternating 
spin chain with antiferromagnetic interactions. These 
two types of ferrimagnets are equivalent to each other in 
certain limits. Many of their thermodynamic or dynamic 
response functions show qualitatively similar behaviors. 
For convenience in the discussion below, we call the fer- 
rimagnets formed from monospins or from alternating 
large and small spins as Type-I or Type-II ferrimagnets. 

The simplest bond-alternating spin chain is a 
ferromagnetic-antiferromagnetic (F-AF) alternat- 
ing Heisenberg spin chain. A typical example of 
this F-AF alternating spin chain is the compound 
(CH3)2CHNH3CuCl3. Manaka et al.^ did birefringence 
measurement on this material. Although the ground 
state of the compound is nondegenerate, a double-peak 
structure of the specific heat was observed. 

A Type-I ferrimagnet can be constructed 
from a tetrameric ferromagnetic-ferromagnetic- 
antiferromagnetic-antiferromagnetic (F-F-AF-AF) 
Heisenberg exchange spin chain. Recently, Hagi- 
wara et al. studied thermodynamic propertieai^ of 



a tetrameric chain compound Cu(3-Clpy)2(N3)2(3- 
Clpy=3-Chloroyridine), abbreviated as CCPA below iS 
Typical behavior of ferrimagnets was revealed. Further- 
more, by exactly diagonalizing an isotropic F-F-AF-AF 
Heisenberg model on small lattice, they found that 
the experimental data for the magnetization and sus- 
ceptibility can be quantitatively understood from this 
modeli^ There were also other theoretical and numeri- 
cal studies on this tetrameric chain system. Yamamoto 
calculated the zero-field specific heat and susceptibility 
of the F-F-AF-AF Heisenberg model using the modified 
spin-wave theory as well as the Quantum Monte Carlo 
(QMC) methodpiS Nakanishi and Yamamotoii found 
that the specific heat of this model shows a double-peak 
structure. However, a direct comparison between 
experiments and numerical calculations, especially for 
the temperature dependence of the specific heat, is still 
absent. 

A typical Type-II ferrimagnet is the alternating S ~ 1 
and 1/2 spin chain. Kahn's groupie synthesized success- 
fully a class of bimetallic chain compounds with each 
unit cell containing two spins with different values. Typ- 
ical compounds include ACu(pba)(H20)3 ■ nH20 and 
ACu(pbaOH)(H20)3 • nH20, where A = Mn, Fe, Co, 
Ni, Zn, pba = l,3-propylenebis(oxamato), and pbaOH — 
2-hydroxy-l, 3-propylenebis(oxamato). This has stimu- 
lated theoretical studies on alternating spin chains. From 
numerical simulations, Pati et al? revealed the coex- 
istence of gapless and gapped excitations in alternat- 
ing spin-(l, 1/2) ferrimagnetic chains using the density- 
matrix renormalization group (DMRG)ii^ The coexis- 
tence of ferromagnetic and antiferromagnetic aspects 
was also found by Brehmer et al—f- from the Quantum 
Monte Carlo (QMC) simulations. Thermodynamic prop- 
erties of the alternating spin chains were studied by the 
transfer-matrix renormalization group (TMRG).^ To fur- 
ther explore the dual features of ferrimagnetism, the 
Schwinger boson and the spin-wave approximations were 
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exploitediiSiia 

In this paper, we investigate numerically thermody- 
namic properties of tetrameric spin chains with alternat- 
ing F-F-AF-AF exchange interactions and compare the 
results to the experimental data for CCPA. In Sec. ^1 
the isotropic F-F-AF-AF Heisenberg exchange model for 
the tetrameric spin chains with its general properties is 
introduced. In Sec. IIIII thermodynamic quantities of the 
tetrameric model including the magnetization, uniform 
and staggered magnetic susceptibilities, and field depen- 
dent specific heat, are evaluated using the transfer-matrix 
renormalization group Section IIVI compares the 
numerical results to the experimental data on CCPA. 
Our numerical results agree well with the experimental 
data for the zero-field uniform susceptibility as well as the 
magnetization. The numerical results for the difference 
of the specific heat between two different fields are also in 
good agreement with the experimental data. We briefly 
discuss the sharp low-temperature peak that appeared 
in the zero- field specific-heat curve published in Refs. 
and0. It is suggested that this experimental sharp peak 
is an extrinsic feature and most probably because of the 
contribution of magnetic impurities. 



II. MODEL HAMILTONIAN AND GENERAL 
PROPERTIES 

CCPA is a tetrameric spin-chain compound. It consists 
of copper chains in which adjacent Cu^^ ions are linked 
by two kinds of azido bridges, namely, end-on and end-to- 
end bridge, respectively [Fig. 1 of Ref. Q- The exchange 
interactions between adjacent Cu^"*" ions are determined 
by the azido bridges linked to them. The low-energy 
magnetic excitations in this material can be modeled by 
a S = 1/2 Heisenberg Hamiltonian defined by 
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where L is the number of unit cells, Jp and Jaf are the 
ferromagnetic and antiferromagnetic exchange constants, 
respectively, g is the magnetic g factor and h is the ex- 
ternal magnetic field. This model belongs to the family 
of Lieb-Mattis-type ferrimagnets, as shown in Fig. ^ 
According to the Lieb-Mattis theorem,^ this tetrameric 
system has {2L + I)-fold degenerate ground states with 
total spin S'tot — L. The macroscopic magnetization 
of the ground state has been quantitatively confirmed 
by magnetization measurements with pulsed and static 
fields on a single crystal of CCPA.^ 

In the limit Jp = 0, the system is decoupled into iso- 
lated trimers, separated by local one-half spins. The 




FIG. 1: Schematic representation of the chain structure of 
Cu(3-Clpy)2(N3)2 (CCPA) with interactions described by 
Hamiltonian Q. The system can be decomposed into two 
sublattices according to the nature of couphng, i.e., A and B 
sublattices, respectively. 



Hamiltonian then becomes 
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This model is exactly soluble. In the absence of an 
applied field, an antiferromagnetic trimer is a three- 
level system with a doubly degenerate ground state of 
S'tot = 1/2. One excited state is twofold degenerate with 
S'tot = 1/2 and an energy Jaf above the ground state. 
The other is fourfold degenerate with Stot = 3/2 and 
an energy 3Jaf/2 above the ground state. For finite 
but small Jp, the ferromagnetic correlation dominates 
the low-temperature excitations. The low temperature 
behavior of the system is expected to act similarly as a 
spin-i ferromagnetic spin chain. 

In the limit \ Jp \^ Jaf, the three ferromagnetically 
coupled neighbors in one unit will bind together to form a 
3/2 spin. In this case, the system is expected to behave 
similarly as a (S, s) = (3/2,1/2) Heisenberg ferrimag- 
netic spin chain. 

For the alternating (S, s) = (1,1/2) Heisenberg spin 
chain there exist both gapless ferromagnetic and 
gapped antiferromagnetic excitations. For the bond- 
alternating ferrimagnetic spin chain considered here, 
these two kinds of excitations are also expected to ex- 
ist. In low temperatures, thermodynamic properties of 
the system are more strongly affected by ferromagnetic 
excitations. However, with increasing temperatures, the 
contribution from gapped antiferromagnetic excitations 
increases. This leads to a crossover in the behavior of 
C and X- Foi' instance, in low temperatures, for the 
(S, s) — (1, 1/2) spin chains, the specific heat and the sus- 
ceptibility vary respectively as C cx T^^^ and x o<; T~'^, 
same as for the spin-^ ferromagnetic spin chain. How- 
ever, in the intermediate temperature regime, they show 
a Schottky-like peak and a minimum, respectively!^ 



III. NUMERICAL RESULTS 

In this section, we will use the TMRG to evaluate the 
temperature dependence of the magnetic susceptibility 
and the specific heat for the tetrameric bond-alternating 
spin model described by Hamiltonian The TMRG 
is a powerful numerical tool for studying thermodynamic 
properties of one-dimensional quantum systems.'iS'iS*^ 
It starts by expressing the partition function of a one- 
dimensional quantum lattice system as a trace of a virtual 
transfer matrix 7j\/ using the Trotter-Suzuki decomposi- 
tion 



Z = Tre-'^^= lim TrT,^/^ 



(6) 



where M is the Trotter number, N is the system size, 
and T = f3/M. For the tetrameric model considered here, 
N — 2L and each site used in the Trotter-Suzuki decom- 
position contains two spins. In the limit N — > oo, one 
can evaluate nearly all thermodynamic quantities by the 
maximum eigenvalue Amax and the corresponding left 
{i/j^ I and right | ip^) eigenvectors of the transfer matrix 
T/\/. For example, the free energy F, internal energy [/, 
and longitudinal magnetization AIz per unit cell can be 
expressed, respectively, as 



F 
U 



— lim — — - In Z 

N~.oo NI3 

(V'^ I fu I V'^) 

Amax 
(V-^ I Tm I ^''^ 



A 



max 



(8) 



(9) 



where the definition of the transfer matrices Tjj and Tm 
can be found from Refs. and ^3- The specific heat 
and magnetic susceptibility can then be calculated by 
numerical derivatives of U and , respectively. 



C 

X 



dU 
dH 



(10) 
(11) 



In our TMRG iterations, 60 states are retained in the 
calculation of the susceptibility, and 80 states are re- 
tained in the calculation of the specific heat. The er- 
ror results from the Trotter-Suzuki decomposition is less 
than 10"'^. The truncation errors are smaller than 10~^. 

Figure|21shows the temperature dependence of the uni- 
form zero-field susceptibility multiplied by temperature 
xT . When Jp = 0, approaches a finite value in the 
limit T 0, in agreement with the Curie- Weiss law for 
a free magnetic moment of S — 1/2. With increasing 
temperature, excitations corresponding to higher spins 
in the trimers are stimulated and increases and satu- 
rates in high temperatures. For finite Jf, diverges in 
the limit T — > 0, because of the formation of ferromag- 
netic long-range order in the ground state. Generally, 




FIG. 2; TMRG results (except for the case Jf = 0) of the 
uniform zero-field susceptibility multiplied by temperature for 
the tetrameric bond-alternating spin model. 



xT decreases monotonically with temperature for ferro- 
magnets, whereas it increases monotonically with tem- 
perature for antiferromagnets. The behavior of xT for 
this ferrimagnet reflects the interplay of these two op- 
posite trends. Gapless ferromagnetic excitations reduce 
the magnetization at low T, while thermally activated 
antiferromagnetic excitations give rise to the increase of 
xT at high T. The presence of a minimum in xT for all 
nonzero Jp is a typical feature of one-dimensional (ID) 
ferrimagnets. Similar behavior was observed in the al- 
ternating (5, s) = (1,1/2) spin chain i4 This minimum 
moves slightly toward higher temperature with increas- 
ing I Jp I / Jaf, resulting from the enhancement of fer- 
romagnetic correlations. 

Figure O shows the temperature dependence of the 
staggered susceptibility. To calculate the staggered sus- 
ceptibility, we add a staggered magnetic field instead 
of a uniform field to the Hamiltonian. Two kinds of stag- 
gered fields are considered. One is a staggered field ap- 
plied to the whole lattice. The corresponding Zeeman 
interaction, staggered magnetization, and susceptibility 
are defined by 
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The staggered susceptibility Xs^'' as defined corresponds 
to the spin structure factor at g = tt. Another kind of 
staggered field we consider is that applied only to the 
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FIG. 3: The zero-field staggered susceptibilities multiplied by 
temperature, xi^' a-nd Xs^' are defined by Eqs. 11411 and (1171 . 
respectively. 



first and third spins in each unit cell. The corresponding 
Zeeman interaction, staggered magnetization, and sus- 
ceptibihty are defined by 



X 



(2) 



4=1 



(15) 
(16) 
(17) 



(2) 

Xs corresponds to a sum of the spin structure factor 
at g = ±7r/2. As revealed by Fig. |2| both staggered 

susceptibihties, Afi^^ and Mi^\ diverge at T = in the 
zero-field limit. 

The divergence of the uniform susceptibility x as well 



as the two staggered susceptibilities xi^'' and Xs^'' is an 
indication of the formation of ferromagnetic as well as 
antiferromagnetic long-range orders at zero temperature. 
This is consistent with a rigorous theorem about the fer- 
romagnetic and antiferromagnetic orders in the ground 
state of Type-II ferrimagnetSflS' 

Figure ^ shows the temperature dependence of the 
zero-field specific heat at different Jp/ Jaf- When Jp — 
0, the specific heat drops quickly in low temperatures be- 
cause of the finite-energy gap between the ground state 
and excitation states. At finite Jp, a small peak devel- 
ops in low temperatures and C shows a small peak-valley 




FIG. 4: TMRG results of the zero-field specific heat at 
I Jf I /Jaf = 0,0.05,0.15,0.5,2.0 in ascending order along 
the direction of arrows. The case | Jf \ /Jaf ~ 2.0 is plot- 
ted as a dashed line for visual effect. The inset shows the 
low-temperature part of the specific heat. 




KbT/Jaf 



FIG. 5: The zero-field specific heat coefficient C/T versus T 
for I Jf I /Jaf = 0, 0.05, 0.15, 0.5, 2.0. 



structure. The low-temperature peak of C results from 
the gapless ferromagnetic excitations, as illustrated by 
Fig. [S] In low temperatures, the temperature depen- 
dence of C/T measures the energy dependence of the 
density of states of low-lying ferromagnetic excitations. 
For finite Jp , since the density of states of ferromagnetic 
excitations diverges at zero energy, C/T is expected to 
diverge at T — OK. For | Jp \ /Jaf ^ 0.5, as shown 
by Fig. 21 the small peak-valley in the low-temperature 
specific heat changes gradually into a hump structure. 
A more conspicuous peak around ksT ~ 0.5Jaf ex- 
ists in all the cases studied. This broadened peak, often 
called the Schottky-like peak, is associated with gapped 
antiferromagnetic excitations, as predicted by various 
methodsSiii^iii for the Lieb-Mattis-type ferrimagnets. 
The double-peak structure of the specific heat for 




ferrimagnets was also observed by Nakanishi and 
Yamamotoii using modified spin-wave theory. They 
found that this double-peak structure of the zero-field 
specific heat is an intrinsic feature with topological origin, 
since the dual features of ferromagnetism and antiferro- 
magnetism in ferrimagnets can potentially induce a low- 
temperature peak as well as an intermediate-temperature 
peak. They calculated the low-lying spectra of the 
tetrameric chain. Three bands of ferromagnetic excita- 
tions with one gapless and two gapped, and one band of 
gapped antiferromagnetic excitation were found. From 
the above qualitative analysis on the F-F-AF-AF model, 
we believe that in the region | Ji? |<C Jaf, the low-energy 
physics of the tetrameric chain is governed by the ferro- 
magnetic coupling. In this region, the peak should move 
to higher temperature with increasing | Jf | • Our numer- 
ical results support the above arguments. The positions 
and the heights of the low-temperature peaks in the cal- 
culations are consistent with the results for the = 1/2 
Heisenberg model with the same ferromagnetic interac- 
tion Jp- With increasing | Jp |, we find that the peak 
moves toward higher temperature and eventually mixes 
with other excitations, leading to a hump structure. 
However, the position and the height of the Schottky- 
like peak are not sensitive to the ratio | Jp \ /Jaf- This 
is different to the modified-spin-wave theory^ It should 
be emphasized that the intrinsic double-peak structure 
of the specific heat is not a distinctive character of bond- 
alternating-type ferrimagnets. As pointed out by Nakan- 
ishi and Yamamoto^ this structure appears when the 
antiferromagnetic gap is larger than the ferromagnetic 
bandwidth. For the (S', s) alternating-spin-type ferrimag- 
nets, this condition is satisfied for S ^ s. 

At finite magnetic field, due to the field-induced split- 
ting of ferromagnetic excitations, a double-peak structure 
of the specific heat can also appear in some ferrimag- 
nets, such as for the case. {S,s) = {1,1/2)M However, 



these two cases should be clearly distinguished. Figure 
shows the temperature dependence of the specific heat 
for Jp/ Jaf — —0.15 at different external fields. When a 
weak field is applied, C shows a sharp low-temperature 
peak. With increasing field, the double-peak structure 
is smeared out and merges together, due to the Zeeman 
splittingiiS The Zeeman splitting of the eigenvalues leads 
to a finite energy gap in the ferromagnetic excitations 
proportional to h, but reduces the energy gap in the an- 
tiferromagnetic excitations. 

Recently, Strecka et alm^^ investigated thermody- 
namic properties of a spin-i Ising-Heisenberg chain 
with F-F-AF-AF bond-alternating interaction using a 
mapping-transformation technique. They considered the 
Ising-type ferromagnetic interaction. Nevertheless, their 
results agree qualitatively with ours. 

IV. COMPARISON TO EXPERIMENTS 

In this part we calculate the susceptibility and specific 
heat of the tetrameric model in the parameter regime 
relevant to CCPA and compare the numerical results to 
the experimental data published in Refs. and|^ 

A. Magnetization and susceptibility 

Figure □ compares the TMRG results of x(T)r for 
different Jf/Jaf to the experimental data for CCPA. 
The parameters used are the same as in Ref. 0- A good 
agreement between experiments and calculations is found 
when Jf/Jaf ~ —0.6, same as in Ref. However, 
above 70 K, the numerical curves are slightly below the 
experimental ones. 

Figure ISl compares the TMRG results of the field de- 
pendence of the longitudinal magnetization to the exper- 
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FIG. 9: TMRG results for the temperature dependence of the 
specific heat ai H = 0, 0.5, 1.0, 1.5, 2.0 T in ascending order 
along the direction of arrows. The parameters used are the 
same as for Fig. |H| The inset shows the low-temperature part 
of C. 



FIG. 8: Comparison of the measured experimental data 
of the field dependence of low-temperature longitudinal 
magnetization^ to the TMRG results. The parameters used 
are Jaf/^s = 29.6 K, Jf/Jaf = -0.6, g = 2.24. is the 
saturation magnetization. 



and field dependence of the magnetization of CCPA with 
Jf/Jaf ~ —0.6. 



imental data. Within experimental errors, the TMRG re- 
sults agree with the experiments in the low-field regime. 
In the high-field regime, the numerical curve deviates 
slightly from the experimental data. It is unknown 
whether this deviation is because of some unknown ef- 
fects or measurement errors. At T = 1.7 K, there is a 
broad magnetization plateau. This plateau is a typical 
feature of a quantum ferrimagnet, due to macroscopic 
magnetization of the ground state. At relatively high 
temperature T = 4.2 K, the plateau region shows a weak 
field dependence resulting from thermal fluctuations. 

At high fields, the magnetization increases sharply 
and becomes saturated forT=1.7K. Atr = 4.2K, 
the magnetization curve is smoothed by thermal fluctu- 
ations. In the work of Strecka et ai^^^ quantum fluc- 
tuations are restricted to the antifcrromagnetic trimers 
since the ferromagnetic coupling is Ising-like. A stepwise 
increase of magnetization with B at very low tempera- 
ture is expected because of the short correlation length 
inherent in this quantum ferrimagnet. Nakanishi and 
Yamamoto calculated the spectrum of the F-F-AF-AF 
Heisenberg model using the modified-spin wavc.^^ They 
found that in the region | Jp \ / Jaf < 1, the antifcrro- 
magnetic excitations are almost dispersionless. It is this 
dispersionless antiferromagnetic excitation that enhances 
the magnetization. The above comparison suggests that 
the F-F-AF-AF Heisenberg model can describe quantita- 
tively the temperature dependence of the susceptibility 



B. Specific heat 

Figure El shows the numerical results on the temper- 
ature dependence of the specific heat C{T) at different 
fields for Jf/Jaf = —0.6. At zero field, C increases in 
a certain power law of T in low temperatures. When an 
external field is applied, a double-peak structure appears. 
With increasing H, the low-temperature peak moves to- 
ward higher temperature, but the high-temperature peak 
moves in the opposite direction. 

Figure [TUI compares the experimental data of C to our 
numerical calculations. At finite fields, the numerical re- 
sults agree well with the experimental data below 3 K. 
However, in relatively higher temperatures, the experi- 
mental curves are clearly above the numerical ones. 

The experimental data shown in Fig. ^|were obtained 
from the measurement raw data by subtracting the 
phonon contribution to the specific heat. This subtrac- 
tion may not always be accurate. Thus we have reason 
to suspect that the deviation at high temperatures may 
come from the error in the subtraction of the phonon 
contribution. 

In order to reduce the error in the data subtraction, 
we compare our numerical results of the difference of 
the specific heat at a finite field and that at zero field, 
ACi = C{T,H) - C(T,0), to the corresponding exper- 
imental data. As shown in Fig. 1111 the agreement 
between the experimental and numerical results is im- 
proved, especially in high T, from 1.3 K to 6 K. How- 
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FIG. 10: Comparison of the specific lieat of CCPA to the 
TMRG results. In (a), the triangle (solid), circle (dashed) 
lines represent the experimental (numerical) data for H — 
0, 0.5 T, respectively; In (b), (c), (d), the circle and solid 
lines represent experimental and numerical results, respec- 
tively. The parameters used are the same as for Fig. |S] 
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FIG. 12: Comparison of the specific heat of CCPA to TMRG 
results after subtracting the corresponding data ai H = 0.5 T. 
Circles and solid lines represent experimental and numerical 
results, respectively. The parameters are the same as for Fig. 
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FIG. 11: Comparison of the specific heat of CCPA to the 
TMRG results after subtracting the corresponding zero-field 
data. Circles and solid lines represent experimental and nu- 
merical results, respectively. The parameters are the same as 
for Fig. H 



ever, in low temperatures, the agreement became worse. 
This is clearly due to the presence of the experimental 
sharp peak around 0.5 K at zero magnetic field. 

Hagiwara et ali suggested that when a Heisen- 
berg spin chain has more than two kinds of exchange 
interactions, the specific heat is expected to show 
more than one peak. Manaka et al& also found the 
double-peak structure in the temperature dependence 
of the specific heat by the birefringence measurements 
in ferromagnetic-dominant F-AF alternating Heisenberg 
chains (CH3)2CHNH3CuCl3. From Sec. we know 



that the intrinsic double-peak structure of the specific 
heat for the Heisenberg F-F-AF-AF model could be ob- 
served when I Jp I / Jaf is small enough. But for 
the given parameter Jf/Jaf — —0.6, the sharp peak 
at 0.5 K is not expected. Strecka et alm^ also made 
a comparison of theoretical results on the F-F-AF-AF 
Ising-Heisenberg model with experiments. In order to 
fit the striking low-temperature peak around 0.5 K, | 
Jf I /Jaf had to be drastically reduced, and the re- 
sults merely qualitatively agreed with experiments. On 
the other hand, from the experimental data, the experi- 
mental sharp peak around 0.5 K seems to be suppressed 
by magnetic field [Fig. 4 in Ref. Q- The behavior is dif- 
ferent from the usual situation, as shown in Fig. El where 
the low-temperature peak of the double-peak structure is 
strengthened by a weak external field. 

We believe that the experimental sharp peak of the 
specific heat around 0.5 K at zero magnetic field is not 
an intrinsic property of the model. It may result from 
some extrinsic properties of the material, such as defects 
or boundary states. Further measurements and more ex- 
tensive theoretical and numerical analysis are needed for 
fully understanding this behavior. 

To avoid the contribution from the experimental sharp 
low-temperature peak of the specific heat at zero field, 
we show in Fig. [HI the results for AC2 = C(T, H) - 
C(T, 0.5 T), which subtracts the data in the field H = 
0.5 T, instead of at zero field. We find that our numerical 
results agree excellently with the experiments. This indi- 
cates that the experimental sharp low-temperature peak 
of the zero-field specific heat is indeed extrinsic. 
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V. CONCLUSION 

We have investigated numerically thermodynamic 
properties of the bond-alternating F-F-AF-AF tetramer 
Heisenberg spin chain using the TMRG. The temperature 
dependence of the spin susceptibility and the specific heat 
is found to be determined by the competition of gapless 
ferromagnetic excitations and gapped antiferromagnetic 
excitations. This leads to a crossover in the behavior 
of C and x- With increasing temperature, xT drops 
sharply in low temperatures, but increases and becomes 
saturated in high temperatures. The minimum of xT 
increases gradually with increasing \ Jf / Jaf\- Both the 
uniform and staggered susceptibilities diverge in the zero 
temperature limit. This suggests that both ferromagnetic 
and antiferromagnetic spin correlations are long-range or- 
dered in the ground state. C shows a small peak in low 
temperatures and a stronger and broadened Schottky- 
like peak in an intermediate temperature regime. The 
nonmonotonic T dependence of xT and the double-peak 
structure of C are also the characteristic features of al- 



ternating spin ferrimagnets. 

Our numerical results for the temperature dependence 
of X and the difference of C between two different 
fields C{H) - C(0.5 T) with Jf/Jaf = -0.6 agree 
excellently with the experimental data for CCPA. The 
field dependence of the magnetization agrees also with 
the experiment. We argue that the sharp peak observed 
at ~ 0.5 K in the zero- field specific heat is an extrinsic 
feature of CCPA. It is likely to be the contribution of 
defects or magnetic impurities. Further experimental 
measurements with high-quality samples are desired to 
clarify this issue. 
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